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Finite size scaling in the solar wind magnetic field 
energy density as seen by WIND 
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Abstract. Statistical properties ol the interplanetary 
magnetic field fluctuations can provide an important insight 
into the solar wind turbulent cascade. Recently, analysis of 
the Probability Density Functions (PDF) of the velocity and 
magnetic field fluctuations has shown that these exhibit non- 
Gaussian properties on small time scales while large scale 
features appear to be uncorrelated. Here we apply the finite 
size scaling technique to explore the scaling of the magnetic 
field energy density fluctuations as seen by WIND. We find a 
single scaling sufficient to collapse the curves over the entire 
investigated range. The rescaled PDF follow a non Gaussian 
distribution with asymptotic behavior well described by the 
Gamma distribution arising from a finite range Levy walk. 
Such mono scaling suggests that a Fokker-Planck approach 
can be applied to study the PDF dynamics. These results 
strongly suggest the existence of a common, nonlinear pro- 
cess on the time scale up to 26 hours. 



1. Introduction 

Statistical properties of the interplanetary magnetic field 
(IMF) fluctuations are a topic of considerable interest in 
space research. The subject is closely related to energy 
transport and acceleration processes in the solar wind [Tu 
and Marsch, 1995; Buriaga, 2001; Goldstein and Roberts, 
1999; Milovanov and Zelenyi, 1998]. An approach is to con- 
sider MHD turbulence as the process responsible for ob- 
served statistical features of the IMF fluctuations such as 
intermittency [Tu and Marsch, 1995; Buriaga, 2001]. Single 
point measurements can not uniquely determine the exis- 
tence of turbulence [Frisch, 1995] and ideally one needs to 
construct a structure function from a range of spatial loca- 
tions in the fluid. However, data taken over long intervals in 
the solar wind is routinely single point and can yield strongly 
suggestive, if not definitive results, and it is this type of data 
that we treat here. 

A hallmark of statistical intermittency is the presence of 
large deviations from the average value on different scales. 
An increased probability of finding such large deviations is 
manifested in the departure of the PDF from Gaussian to- 
ward a leptokurtic distribution [Bohr et al., 1998; Sornette, 
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2000] (see also [van Atta and Park, 1972]). It is well estab- 
lished that such intermittency can be caused by the fluctu- 
ations in the rate of energy transfer of the turbulent cas- 
cade [Neart et ai, 1995] as suggested in Kolmogorov's 1962 
theory [Kolmogorov, 1962]. Following Kolmogorov's ideas 
Castaing proposed an empirical model for the PDF of ve- 
locity differences [Castaing et al, 1990]. Given that: 1) for 
constant energy transfer rate e, the fluctuating quantity has 
a Gaussian distribution and 2) the width of the Gaussian 
has a log-normal distribution it was shown that the resul- 
tant Castaing distribution gives a good fit to the velocity 
difference PDFs in laboratory fluid turbulence experiments. 
In the case of the solar wind, one can obtain a good fit to 
the PDF of velocity and magnetic field fluctuations by this 
method [Sorriso-Valvo et al, 1999, 2001] or can consider the 
fluctuations about the mean value [Padhye et al., 2001], but 
as we shall see this technique can not be applied over the 
entire dynamic range of the magnetic field. Fitting a curve 
to a real PDF over its full range is rather difficult due to 
the large statistical errors present for the less frequently oc- 
curring large fluctuations. In practice one has to discard 
significant amount of data in these tails of the PDFs and 
work only around the center of the distribution curves. For 
example, in Sorriso-Valvo et al. [1999] all events larger then 
3 times the standard deviation of the original sample were 
neglected. The PDF rescaling method used in this paper is 
free of that limitation as we do not fit a specific distribution 
to the data. Instead, we extract the scaling properties of the 
fluctuations directly from the data. Importantly the scaling 
for the entire dynamic range of the investigated PDFs is ob- 
tained from the statistics of the smallest fluctuations. The 
scaling is thus derived using points in the center of the distri- 
bution (with very small statistical errors) and then applied 
to the PDFs as a whole. 

We find that the PDF of fluctuations in magnetic energy 
density is mono scaling. The corresponding collapse of the 
rescaled PDF curves reveals the time scales over which sim- 
ilar physical processes occur and also confirms any correla- 
tions suggested by the inverse power-law form of the power 
spectra. Moreover, it quantifies the asymptotic behavior of 
the distribution of the fluctuations, which is essential for 
constraining turbulence models [Bohr et al, 1998]. Given 
the scaling exponent we can then use a Fokker-Planck ap- 
proach to model the statistical behavior of the time series. 

The WIND solar wind magnetic field [Lepping et al., 1995] 
dataset [Freeman et al, 2000a, b] comprises over 1 million, 
46 second averaged samples from January 1995 to Decem- 
ber 1998 inclusive, calculated from the WIND spacecraft 
key parameter database. The selection criteria for solar 
wind data was given by the component of the spacecraft 
position vector along the Earth-Sun line, X > 0, and the 
vector magnitude, R > 30 RE. The data set includes in- 
tervals of both slow and fast speed streams. Similarly to 
other satellite measurements short gaps in WIND data file 
were present. To minimize the errors caused by such incom- 
plete measurements we omitted any intervals where the gap 
was larger than 15%. The original data were not averaged 
nor detrended. The data is not evenly sampled but there 
are two dominant sampling frequencies: 1/46 Hz and 1/92 
Hz. We use sampling frequency f s of 1/46 as our base and 
treat other temporal resolutions as gaps when the accuracy 
requires it (r < 92 seconds). 
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2. Scaling of the fluctuation PDF 

Differencing techniques date back to the 1920s when 
Richardson conveyed the image of the turbulent cascade and 
proposed the velocity difference as an ideal field for statisti- 
cal studies [Frisch, 1995]. In MHD turbulence the differenc- 
ing method can also be applied to magnetic field magnitude 
measurements. Here, however, we will first investigate the 
scaling properties of the magnetic energy density fluctua- 
tions given by B 2 (t). Let B(t) represent the time series of 
the magnetic field magnitude and B 2 (t) its corresponding 
energy density. A set of time series SB 2 (t,r) is obtained for 
each value of time lag r: 

SB 2 (t,T) = B 2 (t + r)-B 2 (t) (1) 

The PDF P(SB 2 ,t) is then calculated for each time series 
SB 2 (t, r). Figure 1 shows the resultant unsealed PDFs of the 
magnetic field energy density fluctuations for different values 
of the time lag r = 46x2 N seconds where N = 0, 1, 2, 17 
giving a range of time scales between 46 seconds and about 
70 days. We plot the PDF for all values of SB 2 for which 
there are at least 150 samples per bin. 

A generic one parameter rescaling method (the finite size 
scaling approach, see e.g. Somette [2000]; Mantegna and 
Stanley [1995, 2000]) is applied to these PDFs. We first as- 
sume that the distributions can be described by a stable, 
symmetric law (a Gaussian or Levy [Bardou et ai, 2002]) 
within a finite scaling range. We use the peaks of the PDFs 
to obtain the scaling exponent s. This is important as the 
peaks are the most accurate parts of the distributions. Fig- 
ure 2 shows these PDFs rescaled with respect to their cr(r) 
(see e.g. [Castaing et al, 1990; van Atta and Park, 1972]) 
and we see that these also do not collapse onto a single curve. 

Figure 3 shows the peaks P(0, r) of the unsealed PDFs 
(shown in Fig. 1) plotted versus r on log- log axes. We see 
that there is a range of r up to 26 hours for which P(0, r) is 
well described by a power law r~ s with index s = 0.42±0.02. 
This break in the scaling at r ~ 26 hours is consistent with 
that found previously (see e.g Burlaga [2001]; Sorriso-Valvo 
et al. [1999, 2001]). Within this range we now attempt to 
collapse the entire unsealed PDF shown in Fig. 1 onto a 
single master curve using the following change of variables: 

P(SB 2 ,t) = t- s P s (5B 2 t- s ,t) (2) 

If the initial assumption of the stable and symmetric distri- 
butions was correct a single parameter rescaling, given by 
equation (2), for a mono-fractal process, would give a per- 
fect collapse of PDFs on all scales. A self-similar Brownian 
walk is a good example of the process where such collapse 
can be observed (see e.g. Somette [2000]). Figure 4 shows 
the result of the above one parameter rescaling applied to 
the unsealed PDF of IMF energy density fluctuations shown 
in Fig 1. We see that the rescaling procedure (2) using the 
value of the exponent s of the peaks P(0, r), as shown in Fig 

3, gives good collapse of the curves onto a single common 
functional form for the entire range of the data. 

The successful rescaling of the PDFs now reveals some 
key features of the processes governing the magnetic field en- 
ergy density fluctuations. In particular, as we will show be- 
low, one can quantify an asymptotic behavior of the rescaled 
PDFs. A single parameter rescaling strongly suggests that 
the resultant distribution is stable and as such should con- 
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verge to Gaussian or Levy distribution for large fluctuations 
[Sornette, 2000]. The leptokurtic character of the PDFs is 
still present after the rescaling indicating that the Levy dis- 
tribution is more appropriate to model the functional form of 
the curves. For a system where B 2 has an upper, physical or 
instrumental, limit such a distribution is truncated. Impor- 
tantly, truncated Levy distributions are still 1) stable and 
2) go over to Gaussian for sufficiently large time lag r [Sor- 
nette, 2000]. Calculation of the Hurst exponent for the B 2 
time series, by direct examination of the a 2 oc r 2H relation 
as well as the "growth of range" algorithm, confirms such a 
limit does exist - we find a mean exponent H = 0.28 ± 0.02 
up to 26 hours beyond which the variance of the signal sat- 
urates. We then fit a truncated Levy distribution to the 
large SB 2 asymptotes of the rescaled PDFs. This Gamma 
distribution is of the form [see Sornette, 2000]: 

P(8B 2 ) = C±ea;p(-i^!)|JB 2 r< 1+ *> (3) 

Two dashed lines on Fig. 4 show (displaced for clarity) the 
Gamma distribution of the large energy density fluctuations. 
This best fit of (3) to each of the tails of the rescaled PDFs 
is obtained for the power law index a ~ 0.66 and parameters 
C_ = 0.9, C+ = 1 and A = 100. The maximum time scale 
of about 26 hours beyond which the scaling no longer holds 
is in good agreement with previous results [Sorriso-Valvo 
et al, 1999, 2001], (see also Burlaga [2001] and references 
therein). 

It has been found previously [Burlaga, 2001] that the mag- 
netic field magnitude fluctuations are not self-similar but 
rather multi-fractal. For such processes the scaling derived 
from P(0, t) would not be expected to rescale the entire 
PDF as above. To verify this we repeated the rescaling pro- 
cedure for magnetic field magnitude differences SB(t, r) = 
B(t + t) — B(t). Figure 5 shows the result of one parameter 
rescaling applied to the PDFs of the magnetic field magni- 
tude fluctuations. The scaling procedure is satisfactory only 
up to 3 standard deviations of the original sample. This con- 
firms the results of Sorriso-Valvo et al. [1999, 2001] where a 
two parameter Castaing fit to values within 3 standard devi- 
ations of the original sample yields scaling in one parameter 
and weak variation in the other. Attempts to improve the 
collapse by using information in the tails (values \SB\ > 3<j) 
would introduce a significant error in the estimation of the 
scaling exponent s. However, the difference in scaling be- 
tween the PDFs of SB 2 and SB seen here may simply arise 
from the long range correlations present in the magnetic field 
magnitude measurements. This follows since: 

(SB) 2 = B 2 {t + t) - B 2 (t) - 2(B(t)B(t + r) - B 2 {t)) (4) 

The first two terms, B 2 (t + r) — B 2 {t) are the fluctuation 
of the magnetic energy density and as we have seen this 
quantity is mono-scaling. The last term, which includes 
B(t)B(t + t), will however be significant when the autocor- 
relation function of the IMF magnitude time series is non 
negligible. 

The Fokker-Planck equation provides an important link 
between statistical studies and the dynamical approach ex- 
pressed by the Langevin equation [Sornette, 2000]. The 
mono-scaling of the magnetic field energy density PDF, to- 
gether with the finite value of the samples' variance, indi- 
cates that a Fokker-Planck approach can be used to study 
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the dynamics of the unsealed PDF P(SB 2 ,t) in time and 
with respect to the coordinate SB 2 [van Kainpen, 1992]. 
Recently a fractional Fokker-Planck equation has been ob- 
tained that can be used to treat Levy walk PDF dynamics 
[Schertzer et al, 1995]. Alternatively, the anomalous diffu- 
sion of the PDF, consistent with truncated Levy walk, can 
be obtained by introduction of a functional dependence of 
the diffusion coefficient on the "spatial" coordinate (in our 
case SB 2 ). One can then write a Fokker-Planck equation for 
the evolution of the PDF as: 

^-=Ve(P~f(SB 2 )) + VlDP (5) 

where P = P(SB 2 ,t) is a distribution function, 7 is called 
the friction coefficient and D = D(SB 2 ) is a diffusion coef- 
ficient which in this case depends on the magnetic energy 
density fluctuations. The differential operator V/3 = should 
be understood as: 

V ^ = mm < 6 > 

Such a functional dependence of the diffusion coefficient on 
SB 2 leads to the anomalous diffusion resulting in the depar- 
ture from the well known Einstein law relating variance of 
the signal and the time increments Sornette [2000] . 



3. Summary 

In this paper we have applied the generic finite size scal- 
ing method to study the scaling for the magnetic field en- 
ergy density and magnitude fluctuations of the IMF using 
WIND MFI data. This yields the scaling properties of the 
solar wind magnetic field energy density fluctuations that do 
not assume a particular turbulent model. Instead, we use 
the smallest (and best represented statistically) fluctuations 
to obtain a single parameter that gives the scaling exponent 
sufficient to rescale the 46 seconds averaged data over its 
entire dynamic range in SB 2 . The simplest explanation sug- 
gested by these results is that the IMF energy density fluc- 
tuations are governed by a common physical process that is 
self-similar on the temporal scale up to r ~ 26 hours. The 
rescaled PDF follows a non-Gaussian distribution with an 
asymptotic behavior well described by the Gamma distribu- 
tion arising from a finite range Levy walk. Such rescaling 
also indicates that a Fokker-Planck approach can be used 
to study the evolution (in SB 2 space) of the PDF. A simi- 
lar rescaling procedure applied to magnetic field magnitude 
fluctuations confirms more complex, perhaps multi-fractal, 
character of the underlying dynamics in agreement with pre- 
vious results. 
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Figure 1. Unsealed PDFs of the IMF energy density fluc- 
tuations. Time lag r = 2 k x 46s, where k = 0,1,2, ..,17. 
Standard deviation increases with r. Error bars represent 
3a intervals assuming Gaussian distribution of points within 
each bin. 
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Figure 2. PDFs rescaled with respect to their standard 
deviation <t(t). Shown curves correspond to the scaling re- 
gion from figure 3 i.e., r between 46 seconds and 26 hours. 
Rescaled PDFs show no satisfactory collapse. Error bars as 
in figure 1. 
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Figure 3. Scaling of peaks P(0, r) of the IMF energy den- 
sity fluctuations PDFs. Slope of the line is fitted with 95% 
confidence interval. Scaling extends from r = 46 seconds to 
about r = 26 hours. Error bars as in figure 1. 
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Figure 4. One parameter rescaling of the magnetic field 
energy density PDFs. The curves shown correspond to r 
between 46 seconds and 26 hours. Dashed lines represent 
the Gamma distribution of the large fluctuations for the 
PDF of r = 2 4 . Vertical offset was added for clarity. Error 
bars as in figure 1. 
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Figure 5. One parameter rescaling of PDFs of IMF fluc- 
tuations. Shown curves correspond to r between 46 seconds 
and 26 hours. Note the satisfactory collapse to about 3 stan- 
dard deviations of the original sample and poor collapse of 
the tails. Error bars as in figure 1. 



